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Abstract 



The generic Berry phase scenario in which a two-level system is coupled 
to a second system whose dynamical coordinate is slowly-varying is gener- 
alized to allow for stochastic evolution of the slow system. The stochastic 
behavior is produced by coupling the slow system to a heat reservoir which is 
modeled by a bath of harmonic oscillators initially in equilibrium at tempera- 
ture T, and whose spectral density has a bandwidth which is small compared 
to the energy-level spacing of the fast system. The well-known energy-level 
shifts produced by Berry's phase in the fast system, in conjunction with the 
stochastic motion of the slow system, leads to a broadening of the fast system 
energy-levels. In the limit of strong damping and sufficiently low temperature, 
we determine the degree of level-broadening analytically, and show that the 
slow system dynamics satisfies a Langevin equation in which Lorentz-like and 
electric-like forces appear as a consequence of geometrical effects. We also 
determine the average energy level-shift produced in the fast system by this 
mechanism. 
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I. INTRODUCTION 



In the 14 years since its discovery, Berry's phase has proven to be a fruitful development 
in our understanding of the adiabatic limit of quantum mechanics [I]. In his original analysis 
@, Berry considered a quantum system with a discrete, non-degenerate energy spectrum 
whose dynamics is driven by a set of classical parameters R(£) that vary slowly on the 
time scale of the quantum system. He showed that when R(t) was cycled adiabatically 
through a loop in parameter space, and the quantum system was initially prepared in an 
eigenstate \E(0)) of the initial Hamiltonian H(0), the quantum system returned to the 
initial state \E(0)) at the end of the cycle, to within a phase factor exp(i0). This much was 
in agreement with the quantum adiabatic theorem. What was new was that the phase <fi 
contained a contribution <p g whose origin was found to be deeply geometrical ||, and which 
had been discarded in previous treatments of this theorem. <p g is now referred to as Berry's 
phase, and is a functional fl [R(i)] of the loop traced out by R(t) in parameter space. 

Not long after Berry's discovery, it was suggested by Moody et. al. f|] that Berry's phase 
should be observable in nuclear magnetic resonance. It was argued that Berry's phase would 
produce a shift in the energy levels of the resonating spin which would alter the observable 
resonant frequencies. This shift was subsequently observed by Suter et. al. |J. Berry's phase 
has also been observed in a number of other experimental settings [jT|. For our purposes, 
though, the NMR context will be the most relevant. 

The original Berry scenario was generalized in a number of ways in the years following 
his discovery [|l]. For the purposes of this paper, however, the most interesting development 
was the elevation of the slowly varying classical parameters R(t) to the status of dynamical 
quantum variables The total system now consists of a slow subsystem with coordinates 
R(t) coupled to a fast spin-like degree of freedom a. The separation of dynamic time scales 
allows a Born-Oppenheimer treatment of the coupled dynamics. The fast system is again 
found to develop a Berry phase in those states evolving out of an instantaneous energy 
eigenstate due to the motion of the slow system. Berry's phase again produces shifts in 
the fast system energy levels. These shifts are functionals <5i?[R(t)] of the slow system path 
R(t). Geometrical effects also influence the motion of the slow system. They lead to the 
appearance of gauge fields which act on the motion of the slow system [|^-[T0|. The gauge 
fields produce Lorentz-like and electric-like forces in the classical equation of motion of the 
slow system. 

In this paper, we examine a generalization of the Born-Oppenheimer scenario in which 
the slow system is allowed to evolve stochastically. The stochastic behavior is introduced 
by coupling the slow system to a heat resevoir which we represent by a bath of harmonic 
oscillators initially in thermal equilibrium at temperature T. At first, all 3 subsystems will 
be quantum mechanical. Ultimately, we will be interested in the case where the stochastic 
behavior is classical, so eventually we will take the semiclassical limit of the coupled slow 
system/heat resevoir dynamics. The fast system remains quantum mechanical throughout. 
We will examine how the stochastic motion of the slow system influences the Berry phase 
effects discussed above. 

For easy reference, we summarize our principal results. First, because the semiclassical 
motion of the slow system is a random process, the energy-level shifts produced by Berry's 
phase in the fast system's spectrum will cause these levels to broaden. This can be simply 
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understood in the following way. Due to the stochastic force acting on the slow system (see 
below), the boundary conditions imposed on the path R(£) no longer determine a unique 
path. Instead, all possible realizations of the random process Rj(t) must be considered 
which satisfy the boundary conditions. Each realization Rj(t) will produce a level-shift 
(^[R^t)] with a probability P[Rj(t)]. Thus all possible shifts in the original energy E are 
produced, though with differing probabilities. The result is that the original energy level 
E is broadened out by the stochastic evolution of the slow system, in conjunction with the 
Berry phase induced level-shifts. We will calculate the average energy shift 5E produced, 
and the variance associated with the spread of energy shifts about 5E, under appropriate 
restrictions. Second, we show that the semiclassical motion of the slow system is governed 
by a Langevin equation. Geometrical effects are seen to produce the same Lorentz-like and 
electric-like forces in this equation as were found in the original Born-Oppenheimer scenario 
in which the slow system evolved deterministically. We are able to determine the probability 
distribution P[Rj(t)] which we use to calculate the average energy level-shift and broadening 
produced in the fast system. 

The organization of this paper is as follows. We begin in Section |I| by reviewing those 
Berry phase results which are pertinent to the work described in this paper. In Section |II 
we introduce the model system we shall study. In Section [TV] we set-up a path integral rep- 
resentation of its dynamics. We implement the Born-Oppenheimer approximation, watching 
carefully for the appearance of Berry's phase. We also integrate out the unobserved reser- 
voir degrees of freedom. This yields an effective description of the slow system dynamics in 
which only R(t) enters as a dynamical variable. Although the fast system and the resevoir 
no longer appear in this effective description explicitly, they do produce a back action on 
the effective slow dynamics which we obtain while carrying out the above maneuvers. In 
Section [V] we take the semiclassical limit of the effective slow system dynamics. We find 
that the slow system motion is governed by a Langevin equation in which the geometry- 
induced gauge forces appear. We also obtain the probability distribution governing the 
random process R(t). In Section [VI] we use this probability distribution to set-up a gener- 
ating function for the moments of the energy-level shift 5J5[R(t)]. We use the generating 
function to calculate the average energy-level shift and level-broadening, under appropriate 
restrictions. In Section [VI I| , we summarize our results and make closing remarks. Finally, 
in an Appendix, we calculate the spectral density for classical Brownian motion; the goal 
being to determine what restrictions are needed to insure that the stochastic motion of the 
slow system is adiabatic relative to the fast system. 



II. BERRY PHASE PRELIMINARIES 

In this section we collect the Berry phase results which are relevant to the work we shall 
present below. The reader is referred to Ref. |J for a more detailed presentation. 

In the original Berry phase scenario ||, the focus of attention is a quantum system with 
a discrete, non-degenerate energy spectrum. Its Hamiltonian H[H] is assumed to depend 
on a set of classical parameters R which represent an environmental degree of freedom to 
which the quantum system is coupled. The environment is assumed to evolve adiabatically. 
This produces an adiabatic time dependence in the quantum Hamiltonian, H = H[H(t)]. 
The time dependence of the quantum state \ip(t)) is determined by solving Schrodinger's 
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equation using the quantum adiabatic theorem. Towards this end, one introduces the energy 
eigenstates of the instantaneous Hamiltonian H[R(t)}, 



H[R(t)) \E[R(t)) ) = E[R(t)) \E[R(t)) ) 

It is further assumed that the environment is taken adiabatically around a loop in parameter 
space such that R(T) = R(0), and that the quantum system is initially prepared in an 
eigenstate |£'[R(0)]) of the initial Hamiltonian H[R(0)]. The quantum adiabatic theorem 
states that, at time t, the quantum system will be found in the state |i?[R(t)]) to within a 
phase factor, 



\if?(t)) = exp 



ij E (t) - \ f drE[R(r)] 
n Jo 



\E[R(t)}) . (1) 



The second term in the phase of the exponential is known as the dynamical phase and was 
already familiar from previous studies of the quantum adiabatic theorem. The first term 
represents Berry's discovery, and is referred to as Berry's phase, 

= i £ dr (E[R(t))\-^\E[R(t)}) . (2) 

In the cases where Berry's phase is physically relevant, 7^ is non-integrable: it cannot be 
written as a single-valued function of R over all of parameter space. Simon || showed 
that the quantum adiabatic theorem has a line bundle structure inherent in it, and that 
Schrodinger's equation defines a parallel transport of the quantum state around the line 
bundle. Berry's phase is the signature that the associated connection has non-vanishing 
curvature. 

Berry's phase has been observed in a number of physical systems For present pur- 
poses, the nuclear magnetic resonance experiments are the most interesting. Moody et. al. 
Q pointed out that Berry's phase should alter the observed resonant frequencies in NMR. 
In particular, if one were examining the resonance associated with the pair of levels E and 
E', the shift in the resonant frequency Au; would be, 

a 1e(T) -1e>(T) 
A^o = ij 

where T is the period of the oscillating transverse magnetic field H±(t). The resonant 
frequency shift occurs because a shift 5E is produced in each energy-level E by Berry's 
phase, 

SE=*Zp. . (3) 

In the NMR experiments both E and 7^ were time independent. Since 7e(T) is independent 
of the parameterization of R(t), so long as the new parameterization remains adiabatic, 
eqn. will also be correct when 7^ 7^ constant. One simply reparameterizes the time 
t — > if in such a way that TE(i') is time independent. 

The above scenario can be generalized. We promote R{t) from a classical degree of 
freedom with no dynamics of its own, to a fully-dynamical quantum variable |||7[]. As 
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above, R(t) couples to a spin-like degree of freedom er. For the remainder of this paper we 
will assume that er corresponds to a pseudo-spin 1/2. To stay within the context of the 
quantum adiabatic theorem, we assume that a Born-Oppenheimer treatment of the coupled 
dynamics is appropriate, and that R (er) is the dynamical variable of the slow (fast) system. 
The Hamiltonian for the coupled system is taken to be 

H tot = ^ + V[R]-gR-cr . (4) 

Here P is the momentum conjugate to R; crj are the Pauli matrices; we allow for the 
possibility of a potential V acting on the slow system; and g is a coupling constant. 

In applying the Born-Oppenheimer approximation, one first considers the fast system, 
treating the slowly-varying R as fixed. The fast motion is governed by Hf = —gH ■ er with 
eigenstates |.E±(R)) and energies ±gR, where R = |R|. In fact, R varies adiabatically. From 
the quantum adiabatic theorem we know that when the fast system is prepared initially in 
the state |i?[R(0)]), its state at time t is given by eqns. ([[]) and (0), and Berry's phase 
produces the energy level shift SE given by eqn. @. The slow dynamics is governed by 



H eff = (E(R)\H tot \E(R)) 
(P-A[R]) 



\2 



$(R) + £(R) . (5) 



2M 
Here, 

A(R) = ift(£(R)|V R |£(R)> , (6) 



and, 



2 



$(R) = - $>i(R) . (7) 

z i=± 

From eqns. @ and @, we see that A[R] is related to Berry's phase: hjE = R • A[R]. 
c/y(R) is the quantum metric tensor PJTl| and corresponds to the real part of the quantum 



geometric tensor 2^-, 

^• = (^1(1-1^)^1)1^) 

where 9j = d/dRi. We see that geometrical effects produce gauge potentials $ and A in 
the effective Hamiltonian that governs the slow dynamics. In the semiclassical limit, the 
equation of motion for the slow system is 

MR = R x B[R] - V R [$[R] + V[R] + E[TL]] , (8) 

where B = V x A. Thus, geometrical effects lead to the appearance of Lorentz-like and 
electric-like forces that act on the slow system in the semiclassical limit. 

This concludes our summary of the Berry phase physics we shall need below. 
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III. THE MODEL 



As mentioned in the Introduction, we would like to explore how the Berry phase physics 
discussed in Section |I| is affected by stochastic motion of the slow system. To produce 
stochastic behavior, we introduce a heat resevoir that couples to the slow system, and which 
initially is in thermal equilibrium at temperature T. Because the resevoir degrees of freedom 
are unobserved, they must be traced out of the dynamical description. The result of this 
operation is that stochastic and frictional forces appear in the dynamics of the slow system. 

To produce a tractable model, some restrictions must be placed on the heat resevoir, and 
on its coupling to the slow system. These restrictions, and the model to which they lead, 



have been discussed in great detail in the literature ||12| , [13| . Consequently, our discussion 
here will be brief. The reader is referred to these papers for further discussion. 

The first restriction imposed is that the heat resevoir contain an infinite number of 
degrees of freedom. The idea here is that any energy transferred from the slow system to 
the infinitely many resevoir degrees of freedom effectively disappears into the resevoir, not 
to return in any physically reasonable amount of time. The result is dissipation of energy 
and irreversibility in the motion of the slow system. 

The second restriction is that the slow system couple weakly to the individual resevoir 
degrees of freedom. As such, each resevoir degree of freedom is only weakly disturbed from 
equilibrium. This restriction allows us to analyze the excitation of the resevoir away from 
equilibrium using a harmonic approximation [|T!2| , P~3"| 1 ■ Consequently, the resevoir degrees of 
freedom can be modeled by harmonic oscillators that represent the normal modes of the 
resevoir. 

This combination of restrictions is known to produce stochastic and frictional forces in the 
semiclassical motion of the slow system once the resevoir has been traced out |r^Jl5[| . Because 



Berry's phase requires the slow system to evolve adiabatically relative to the fast system, we 
must add a further restriction to insure adiabaticity of the stochastic motion. The essential 
point is the following. Because of the stochastic force, R(t) contains fluctuations or noise. 
The range of frequencies Auj present in the fluctuations can be determined by examining the 
spectral density Jr(cj) of R(t). In the Appendix we show that, in the semiclassical limit, 
Jr(c<j) is proportional to the spectral density Jf(u) of the stochastic force F(t). The range 
of frequencies (bandwidth) present in F(t) is of order l/r c , where r c is the correlation time 
of F(t). The correlation time is known for our model of the resevoir 



Thus the bandwidth Alu of the random process R(t) is also of order l/r c . We can insure 
that R(t) evolves adiabatically if we require that hAu be much less than the energy-level 
spacing of the fast system, AE = E + — E- . This is the final restriction we need to impose 
on the slow system/reservoir dynamics. For example, since we see that when 

IK < T < 100 K, the bandwidth for fluctuations in R(i) satisfies 10" 4 eV < hAcu < 10~ 2 eV. 
Thus if AE > 0.1 eV, the slow system will evolve adiabatically relative to the fast system, 
even though its dynamics is stochastic. 

Assuming that all 3 of these restrictions are satisfied, the following model Hamiltonian 
will produce adiabatic stochastic motion of the slow system: 



H= — + ^[ R ] -0R-o- 
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+ Eqr-q, + r 2 eA • (9) 

j j -'"j^j 

The first term on the right-hand side (rhs) is the non-interacting slow system Hamiltonian 
H s . The second term on the rhs is the slow system/fast system interaction Hamiltonian 
H s f. The third term is the non- interacting resevoir Hamiltonian H r , and the last two terms 
give the slow system/resevoir interaction Hamiltonian H sr . mj and ujj are the oscillator 
masses and frequencies, respectively; and Cj are the coupling constants for the slow sys- 
tem/oscillator interactions. These constants are not determined by the model, though they 



satisfy a constraint |12| that connects them to the friction coefficient rj that will appear in 
the Langevin equation obtained in Section |V[ The last term on the rhs is a counter-term 
introduced to insure that ^(R) is the potential acting on the slow system in the semiclassical 
limit. 

Having introduced our model, our next task is to set-up a path integral treatment of its 
dynamics. 

IV. DYNAMICS VIA PATH INTEGRALS 



In this Section we apply the influence functional formalism of Feynman and Vernon [13 
to describe the dynamics of our model. The analysis begins with the density matrix for the 
composite system (slow system/fast system/resevoir). The essentials of its time development 
are presented in Section PV A[ The Born-Oppenheimer approximation is implemented in 
Section [IV B| , and the resevoir degrees of freedom are traced out in Section |1V l| . The main 
result of this section will be a path integral representation of the reduced density matrix 
describing the effective dynamics of the slow system. Geometrical contributions to this 
effective dynamics appear during implementation of the Born-Oppenheimer approximation; 
stochastic and dissipative effects enter in when tracing out the resevoir, though this will not 
become apparent until Section |V|. 



A. Density Matrix Preliminaries 

We begin with the density operator p(t) of the composite system. Its evolution is given 



by, 



pit) = exp 



/> ( -— j X (>xp 



IK' 1 



(10) 



H is the Hamiltonian for the composite system (see eqn. and the time evolution begins 
at t = -T /2. 

It proves useful to work in the basis |x, a, z ) in which the slow system is at x, the fast 
system has pseudo-spin projection o along x, and the harmonic oscillators have positions 
Z = (. . . , Zj, . . .), where j labels the oscillator degrees of freedom. In this basis, eqn. (p~0|) 
becomes, 
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{xf,a f ,Qf\p(t) \y f ,a' f ,Z f ) = 

J dx'dy'dQ'dZ' (x>,Q'|p(-T/2)|yV,Z'> 

x K(pc f , 07, Q/, t- x', a, Q', -T/2) iT(y /; a), Z,, t; y', a', Z', -T2). (11) 

Here, 



if(x/, or/, Q/, t; x', tr, Q', -T/2) = (x /; o-/, Q/| exp 



|x',<t,Q'>. 



if* is the complex conjugate of if and will not be written out explicitly. In Section [IV B| we 
will obtain a path integral expression for K. 

We must specify an initial condition for the time development. We assume that the 
resevoir is uncoupled from the slow/fast systems prior to t — —T/2. Thus p(— T/2) factors: 



p(-T/2) = p r (-T/2) p, f (-T/2) 



(12) 



The resevoir is assumed to be in thermal equilibrium at temperature T initially. Thus the 
energy-levels of each oscillator degree of freedom are populated according to the Boltzmann 
distribution, and p r .(— T/2) is the product of the density matrices of the individual oscillators 
T3fl . The initial condition on p s /(— T/2) is that it correspond to a Born-Oppenheimer pure 



state. Thus the slow system is prepared in the position eigenstate |x'), and the fast system 
in the negative energy state l-ELfx']) corresponding to the pseudo-spin aligned along x' (i. e. 
a = 1/2): p a f{-T/2) = |x', l/2)(x', 1/2| = p E _{-T/2). Thus, 

{x f ,a f ,Q f \p(t)\y f ,a' f ,Zf) = 

J dx'dy'dQ'dZ' p r (Q', Z\ -T/2) ~p E _ (x', y', -T/2) 

x K(x f , <r f , Q f , t; x', 1/2, Q', -T/2) K*(y f , a' f ,z f , t; y', 1/2, Z', -T/2). (13) 



B. Born-Oppenheimer Treatment of Slow/Fast Dynamics 

Implementing the Born-Oppenheimer approximation in Berry phase systems is well- 
understood [|16| -|T8lj. We present the analysis for K. A similar analysis applies for K*, but 
will not be given here. 

K describes evolution over the time interval [—T/2, t}. We introduce intermediate times 
ti, (i = 1, . . . , N — 1), and break the evolution into N smaller intervals of duration e = 
(t + T/2)/N. Thus, 

N 

K = ( X/ , a f , Q/| I] U(U, i.-OIx', 1/2, Q'), (14) 

with U{ti,ti_i) = 1 — ieH(ti)/h. Inserting basis sets |xj,o~i,Q-) = |xj) l-E^Jxj] )|Q-) at the 
intermediate times ti, eqn. (0) breaks up into a product of factors describing evolution over 
the time intervals (tj-i, U). The factor associated with the z-th time interval is, 



S 



(xi,(Xi,Q-| 
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1 — -jr(H s + H s f + H r + Hg 



I x i-l! Ci-l, Qi-l) 



(xi.Qil 



l^Qli) (15) 



Recall that the slow system evolves adiabatically relative to the fast system. Thus the fast 
system will remain in the initial (negative) energy-level (i. e. a,; = 1/2 for all i). Furthermore 



and, 



^.[Xip^jx^] ) = 1 + iey-fe) + 0(e 2 ), 



Kjx 2 ]| H s \E ai _^ t ] ) = H eff (U) + 0(e 2 



Here 7-(t) is the Berry phase associated with the fast system state |-E-[x] ), and H e ff(t) is 
given by eqn. (|5|), with the state appearing in eqns. @ and (0) given by |-E-[x] ). 
Inserting all these results into eqn. ([14]) gives, 



K = ( x /7Q/|exp 



fl J-T/2 



drH'(r) 



|x',Q'), 



(16) 



where, 



T 



-fry_(r) + E.(t) + H eff (r) 



Notice that explicit reference to the fast system has disappeared because the adiabatic time 
dependence and the initial condition force it to remain in the £L energy-level at all times. 
In its place, the back action of the fast system on the slow system has produced the Born- 
Oppenheimer potential energy -E_[x], and the gauge potentials $[x] and A[x] in H'{r). 
Equation ( |16D can be written as a path integral in the usual way ||19|| , 



K(xf,(T f ,Q f ,t;x!,l/2,Q , ,-T/2) = 6 l , fA/2 J Px(t)X>Q(t) exp 
Here, 



t(.^s + S r + S s 
a 



S s = Jdt —± 2 + x • A[x] - $[x] - £_[x] - V[x\ 

s r = /^£^(Q?-^Q 



2 ) 



dt 



ec,x,.q,-x 2 x: 



. 2 mj u] 



(17) 

(18) 
(19) 

(20) 



All paths x(t) (Q(i)) appearing in eqn. ( ]T7| ) begin at x' (Q') and end at xj (Q/). 
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C. Tracing Out the Resevoir 



When K and K* are substituted into eqn. flT3|), we find, 
(x f ,(jf,Q f \p(t)\y f ,(j' f ,Z f ) = 8 aftl/2 

x J dx'dy'p E _(x',y',-T/2) 
x | dQ'dZ'p r (Q', Z', -T/2) 

x Jvx{t)Vx'{t)VQ{t)VQ'{t) exp | (S[x, Q] - S[x',Q']) . (21) 

Since the resevoir degrees of freedom are unobserved, we must trace over them. Setting 
Q/ = Zj in eqn. (|2"ID and integrating over Qj, the left-hand side (lhs) becomes the reduced 
density matrix p E _(/Xf,yj,t), while the full equation determines its time dependence, 



pE_(x/,y/,t) = J dx.'dy'J E _(x f ,y f ,t;x.',y', -T/2) p B _(x',y', -T/2), 



where, 



J E _(x f ,y f ,t;x!,y',-T/2) = J Z>x(t)Px'(t) exp 



l(5.[x]-5.M) 



Tjx,x'l. 



(22) 



(23) 



All paths x(i) (x'(t)) appearing in eqn. (|23|) begin at x' (y') and end at xj (y/). T[x, x'] is 
the influence functional which contains all the effects of the resevoir on the motion of 
the slow system. For a resevoir composed of harmonic oscillators, T[x, x'] can be evaluated 
exactly [p~3|j , 



T[x, x'] = exp 



1 rt 

— — / drds [x(r) — x (r) [a(r — s)x(s) — a*(r — s)x'(s)] 
n J -t/2 



tl J -T/2 

and a(r — s) = a R {r — s) + ia/(r — s) (see below). Thus 
J E _(xf,yf,t;x',y',-T/2) = J Vx(t)Vx! (t) 

x exp 



(24) 



5 s [x] - 5 s [x'] - y rfrrfs [x(r) - x'(r)] aj(r - s) [x(s) + x'(s)] 
1 ft 

- / drds [x(r) - x (r)] a fl (r - s) [x(s) - x (s)] 
ft Jr/2 



x exp 



where, 



oiAt — s) 



C 2 

E ~ — " — c °th 



Tiuij " 
2AT 



cos [a^r - s)] 



-E 



sin [vj(T-s)\ 



(25) 

(26) 
(27) 



Eqns. (25)-(p7|) give the effective dynamics of the slow system. Geometrical effects have 
produced gauge fields in S s \pc]] and the stochastic and dissipative effects due to the resevoir 
lurk in the terms containing and aj, respectively, though this will not become apparent 
until we take the semiclassical limit of the effective dynamics. 
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V. SLOW SYSTEM EFFECTIVE DYNAMICS: SEMICLASSICAL LIMIT 



Within the context of our model, we have established the slow system's effective quan- 
tum dynamics. In this Section we will see that the semi-classical limit of this dynamics is 
dissipative and stochastic. This limit is taken in Section |V A| . In Section [VB| we obtain the 
Langevin equation governing the stochastic motion, and the probability distribution func- 
tional characterizing the statistical properties of this motion. The Langevin equation is found 
to contain the same geometric forces that appeared in the deterministic Born-Oppenheimer 
scenario discussed in Section [II]. 



A. Semiclassical Limit of Je_ 

Eqn. (^) describes the slow system effective dynamics. In the absence of a resevoir, the 
influence functional JF[x, x'] = 1, and the semiclassical limit of the remaining exponential 
factor can be obtained by the method of steepest descent. When the resevoir is present, 
jF[x, x'] 7^ 1, and we must examine the form it takes when h — > 0. For a resevoir composed 
of harmonic oscillators, this limit has been carried out in Ref. 



TM. The result is, 



JEAxf,yf,t;x?y,-T/2) = J v*{t)v*!{t) 

S a [x) - S s [x'] - | f ^ dr [x(r) - x'(r)] ■ [x(r) + x'(r)] 
' ,k T f rfr{x(r)-x'(r)} 5 

J-T 2 



x exp 



x exp 



h 2 



(28) 



The physical significance of r] will become clear in Section |V B . 

The real exponential factor in eqn. (p8f) introduces important simplifications when h — > 0. 
In this limit, it is very sharply peaked about x(r) = x'(r) so that the dominant contribution 
to the integral comes from pairs of paths for which |x(r) — x'(r)| <C 1, for all r. Because the 
paths x(r) and x'(r) are continuous, JJf_ is also sharply peaked about x^ = y^ and x' = y'. 
By eqn. (|2~2"D, Pe_(x, y> i s similarly sharply peaked about x = y. Thus, expanding the 



argument of the complex exponential in eqn. ([28 ) to second order in [x(r) — x'(r)] introduces 
negligible error when h — > ||15|| . To carry out the expansion, it proves convenient to 
introduce center-of-mass and relative coordinates, respectively, R(r) and p(t): 



R( 



x(r) +x'(r) 



T 



p(r) = x(r) - x'(r) 



Using eqn. (IB), it is straightforward to show that 



-T/2 



p{r) ■ (-MR + R x B[R] + E[R]) + O (|p 



(29) 



(30) 



Here B[R] = V x A[R]; E[R] = -V ($[R] + V[R] + ^_[R]); and A[R] and $[R] are the 
geometric gauge potentials discussed in Section ||. It follows immediately from eqn. ( |29|) 
that, 
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7/ 



rfr [x(r) -x'(r)] ■ [x(r) +x'(r)] = -77/ drp(r) -R(r). 

2 J-T/2 J-T/2 

Putting together all these results gives, 

J S E C _ (x /)y/ ,t;x',y',-T/2) = j VR(t)Vp(t) 



(31) 



x exp 



/i J-T/2 

x exp 



drp(r) 



-MR — ?,R + R x B + E 



r/2 



dr [p(r) 



(32) 



B. Appearance of Stochastic Dynamics 

To bring out the stochastic character of the dynamics contained in eqn. (|32]), we again 
focus on the real exponential in this equation. Its argument can be re-written as, 



7]kT n 
h 2 J -r/2 



drp(r) ■ p(r) = -—2 / / drds pii^Aijir - s)pj(s), 
-r/2 2n J -r/2 j-r/2 



(33) 



where Ajj(r — s) = 2i]kT5ij5(T — s). To bring out the meaning of this term we introduce a 
Hubbard-Stratonovitch transformation: 



$[p] = exp 
1 



2h" J -r/2 ./-r/2 
£>F(t) exp 



-2 / / dTdspi{T)Aij{T - s)pj(s 
ft J-T/2J-T/2 

\LnJL dTdaF - iT)A * {T - s)F > ia) 



x exp 



J-r/2 



drp(i 



F(r) 



(34) 



where ^(r-s) = (l/2r]kT)5 ij 5(r-s), and JV is an (infinite) normalization constant which 
we suppress below. Eqn. (0) is an identity which can be proved by evaluating the Gaussian 
integral on the rhs. From eqn. (|3~4]), it is clear that $[p] is the characteristic functional 
for the Gaussian random process F(t). Functional derivatives of $[p] with respect to p(t) 
generate the correlation functions of F(t). In particular, 

2 



<^W-(s)> 



h 



5 2 $[p] 



ij 5pi{r)5pj{s) 



Aj(r - s) 



2r]kT5ij5(r — s) 



(35) 



Thus F(t) has a two-point correlation function identical to that of a classical stochastic 
force produced by a heat resevoir with friction coefficient i] and temperature T. We shall 
see below that this is the correct interpretation for F(t). 
Making use of eqn. 



in eqn. (|32|) gives, 
J s E c _(x f ,y f ,t;x',y',-T/2) = J VR(t)Vp(t)VF(t) 

x exp -— — f dT[F(r)] 2 



exp 



h J-T/2 



drp(r) ■ LR(r) 



(36) 
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where LR(r) = —MR — rjR + R x B + E. We recognize the p(t) -integral as the Dirac delta 
functional. Thus, 



J%_(xf,y f ,t',*,y',-T/2) = J VR(t)VF(t) 



x exp 



Ar/kT J-t/2 



dr[F(r)Y 



LR(t) + F(t) . (37) 



Clearly, in the semiclassical limit, the only paths that contribute to JJf_ are those which 
satisfy the Langevin equation, 



or, 



LR(t) + F(t) = 0, 



MR = -r/R + R x B[R] + E[R] + F(t). 



(38) 



(39) 



As promised in Section ^TTJ, introducing the resevoir has caused the semiclassical dynamics 
of the slow system to become stochastic, rj can be interpreted as a friction coefficient, and 
F(t) as a Gaussian stochastic force. We also see that the geometric Lorentz-like and electric- 
like forces which appear in the deterministic Born-Oppenheimer scenario (see eqn. (§)) also 
appear in the (adiabatic) stochastic generalization of this scenario. Eqn. (|39|) corresponds 
to one of the two main results of this Section. 

To obtain the second, we carry out the F(t)-integration, 



J%_(B. f ,p f = 0,t;B!,pf = 0,-T/2) 



VR(t) exp 



Ar/kT J-T/2 



dr 



LR{t) 



• (40) 



Inserting this into eqn. (|22[) , we find 

p£.(R/,P/ = (M) = 

dR' [ VR(t) exp 



1 



ArjkT J-t/2 



dr 



LR{t) 



p^_(R,p' = 0,-T/2). (41) 



We see that p s £_(Rf,t) is found by summing over all paths which lead from R' to R/, and 
then integrating over all possible R'. From eqn. (f4T|) , the probability that a path R(t) lies 
in a "volume" T>R(t) in the space of paths that join R' to Rf is clearly, 



VP = exp 







2 


f- 1 /' * 


LR(t) 




Ar/kT J-t/2 





VR{t). 



(42) 



This identifies the probability distribution functional P[R(t)] as, 

P[R{t)\ = exp 







2 


f- 1 /' * 


LR(t) 




Ar/kT J-t/2 







(43) 



We will use eqn. (|43|) in the following section to calculate the average shift and broadening 
of the fast system energy-levels. 
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VI. FAST SYSTEM ENERGY LEVEL-SHIFT AND BROADENING 



In this Section we calculate approximately the average shift and broadening produced in 
the fast system energy levels by the combination of Berry phase and stochastic effects. In 
Section [VI A| , we introduce the generating function that is the basis of our calculation, and 
state the conditions we must impose to make the calculation tractable. The generating func- 
tion is evaluated in Section VI B . In Section [VI Q and |VID , we determine the approximate 
average level-shift and level-broadening, respectively. 



A. Preliminaries 

We restrict the slow system to two-spatial dimensions: R(t) = R(t) [cos 4>{t) , sin0(t),O]. 
H s f takes the form: 

/ o e -i<i> \ 

H sf = gR{t) lit o J ' 

We assume throughout this Section that the fast system is initially prepared in the state 
|£-[Ri] )> where R (- r / 2 ) = R Usm § e q ns - © and (0). 

$ _ [R] = (5Zf = *L . A[R] = = R R x 4 . (44) 

1 J 2M 2MR 2 ' 1 J 2 2R 2 K ' 

From these, B = V x A, and E = — V($_[R] + V[R] + £L[R]). To simplify the following 
calculation, we choose V[R] = — £L.[R]. 
Our Langevin equation becomes, 

MR = -77R - ^<5 2 (R)R x z + ^R + F(f). (45) 

The Lorentz force is seen to vanish everywhere except at the origin, while the electric force 
is divergent there and repulsive. Consequently, any finite energy motion of the slow system 
must have a turning point with R tp 7^ 0. In 2D, then, the Lorentz force disappears from 
eqn. (|4"5|). The electric force is seen to cut off rather quickly, and is second order in Planck's 
constant h. Thus, away from the origin, the electric force is small in the semiclassical 
limit. To simplify the following analysis, we replace the exact electric force with something 
qualitatively similar. To capture the strong repulsion at the origin, we simply remove the 
origin from the xy-plane so that R(t) 7^ for all t. To mimic the rapid cutoff away from 
the origin, we simply drop the electric force from eqn. (JR]). With this simplification, our 
Langevin equation reduces to free Brownian motion in the punctured xy-plane: 

LR = MR + 77R = F. (46) 

The boundary conditions are: R(— T/2) = R. ; , and R(T/2) = Rj. 

We write R(t) = Ro(£) + y(t). Here Ro(t) is the homogeneous solution of eqn. ([46"D; 
and y(t) is a particular solution. Ro(£) satisfies the same boundary conditions as R(t), so 
consequently, y(±T/2) = 0. We want a condition that insures that the fluctuating particular 
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solution y(t) remains small compared to the noiseless homogeneous solution Ro(i) in an 
average sense: (|y(t)| 2 ) <C |Ro(t)| 2 . For Brownian motion, it is well-known that ||20|| . 



(lyWI 2 ) 



2kT 
i 

V 



(47) 



Defining R = (Rj + R/)/2; requiring (|y(£)| 2 ) <C R ; and using eqn. P7|) gives the condition 

we seek, 



r]R 2 

r ' 



(48) 



Brownian motion which satisfies eqn. ( f48|) will be referred to as low-noise Brownian motion. 

The central object of this Section is the generating function F(p) for the moments of the 
energy-level shifts 5i?_[R] (see eqn. @): 



(49) 



F{p) = I VR(t) P[R] exp[p5E[R]]. 
It follows immediately that, 



d[\nF(p)} 



dp 



p=0 



JVR(t)P[R]5E[R] w 
JVR(t)P[R] ~ ' 



(50) 



and, 



^™ U = (^-(JI;)W. 



(51) 



We will use the standard deviation a as a measure of the energy-level broadening. 



Using eqns. 



B. Generating Function 



@, ©, @, and (H) gives 



F(p) = / PR(t) exp 



1 / >r / 2 r - • 



(52) 



where e = —2r]kTh/T. It is not possible to evaluate eqn. (]5^) analytically under arbitrary 
conditions. The calculation is analytically tractable if we restrict ourselves to low-noise 
Brownian motion. As discussed in Section |V1 A| , in the low- noise limit, the fluctuating 
component y(t) is always small compared to the noiseless component Ro(t).This allows us 
to expand the argument A of the exponential in eqn. ( ^|) to second-order in y(£), yielding a 
Gaussian path integral. The presence of R 2 in the denominator of the second term appearing 
in the integrand of A complicates the analysis however. Since the noise component y(t) is 
small, R 2 ~ R throughout the motion. Consequently, we will approximate R 2 by R in this 
denominator. Strictly speaking, we are discarding terms that should be kept by doing this. 
As such, our results in Section |VI C| and |VI D| should be thought of more as rough estimates, 
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than as rigorous results (see the discussion at the end of Section [VI D|) . In principle, one 
could repeat the calculation below including these extra terms, though we will not do so 
here. Carrying out this expansion gives, 



Here, 



F(p) = exp [p5E4R Q ]] J Vy(t) exp 



J[y]= M 2 (y) 2 + r ? 2 (y) 2 + 2Mr ? y-y 



1 



ArjkT J-t/2 



r/2 



drl[y) 



epe 3lj 

+ -=2- 

R 



2R o Vj + Vi Vj 



(53) 



(54) 



and e^ k is the Levi-Civita density. 

It proves useful to Fourier transform y(t). A sine-transform is needed since y(±T/2) = 0, 



n b 



y(*) = Ey n sin (?f0 ■ 

n=l V 1 ' 

The need for adiabatic noise requires an upper cutoff Ub = libit /T on the noise spectrum 
(see Section p| ). The boundary conditions produce a low frequency cutoff u c = ir/T. We 
assume T is large, though finite, so that y(t) can be written as a complex Fourier integral, 



y(t) 



2tt 



dujyiuj) exp[— iut], 



(55) 



where continuation to negative a; requires y(— c<j) = —y(u). Note that this is consistent 
with the boundary condition requirement that there be no zero-frequency mode in the noise 
spectrum. Substituting eqn. ( |55| ) into eqn. (|53|), and completing the square in the integrand 
gives (eventually), 



F(p) = exp 



x J Vz(lu) exp 



duj 



ArjkT J-Ub 2-7T 



(56) 



where, 



Carrying out the Gaussian integral gives, 



Y( P ) 

x(pV 



(57) 
(58) 



(59) 



where, 
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and, 



Cl 

c 2 
c 3 



X(p) 



5E4R ] 

T]kTh 2 /T 2 T? 

MR 2 /2rfkT 



exp 



cip + 



c 2 p 2 r^b du 2 f A - p/2 



^6 27T 1 / 3 - p 2 f 4 



/ri c 3(/3-P 2 /4), 



(60) 
(61) 



f 1 = \R {u)\ 2 {M*<J* + 7i 2 u 2 ) 

/a = iw [z • R (cu) x R*(cj)](r/A;T^/r J R 2 ) 

/ 3 = (M 2 c^ 4 + rfu 2 ) 2 

ft = {r]kThuj/TR 2 ) 2 . 



(62) 



C. Average Energy-Level Shift 



From eqn. (|50j 



^- = >np)iu-|;[in^)]u 

= ci - 
= 5£_[Ro]. 



(63) 



In the low-noise limit we find that the average energy-level shift is given by the shift produced 
by the noiseless component of R(t). The absence of a further noise correction to eqn. ((63|) 
is surely a consequence of the low-noise approximation. Such a noise correction is expected 
to occur for stronger noise. 



D. Average Energy-Level Broadening 



From eqn. 



a 



C-2 



IU C IT h JUJ C TT f 3 



(64) 



Solving the homogeneous Langevin equation for Ro(£), and Fourier transforming to ob- 
tain R (c<j); introducing a dimensionless frequency x = u(M/t))\ and using eqn. ( |6~2"D gives 
(eventually), 



a 



h 


(MkT\ 


"IR/-R,! 2 kT 


f x b dx 1 


n 


V v 2 R 2 ) 


[ R 2 ' (r/if/T). 


Jxc TT X 2 (l + X 2 ) 2 



(65) 



A simpler relation can be found if we choose boundary conditions such that, 
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R/ - Rj r, 

which corresponds to fixing the separation of boundary points to be roughly the same size 
as the diffusion cloud for Brownian motion over a time T . Eqn. fl65|) becomes, 




a ~ kT ( ) K ) (66) 

where, 




b dx 



7T X 2 (l + X 2 ) 2 



and r = M/77. As mentioned above, eqns. ( |65D and fl66|) should be thought of more as 
rough estimates, than as rigorous results. Clearly though, we do see that level-broadening 
is produced by the combination of Berry phase and stochastic effects. A Monte Carlo 
evaluation of the path integrals appearing in eqns. ( |50|) and (pl|) would be very interesting 
since such a calculation would not require the simplifications we found it necessary to make 
to produce an analytically tractable calculation. It is important to keep in mind that such 
a Monte Carlo calculation must still satisfy the third condition imposed in Section [ITJ to 
insure that the noise respects the adiabatic requirements of Berry's phase. 



VII. CLOSING REMARKS 

In the usual Berry phase scenario one considers a pair of interacting systems with 
vastly different dynamical time scales, and treats the coupled dynamics using the Born- 
Oppenheimer approximation. In this paper we generalize this scenario, allowing the slow 
system dynamics to be stochastic as well as adiabatic. We introduce a model that allows us 
to study how the usual Berry phase scenario is modified by the stochastic dynamics. 

Our principal results are: (1) a broadening and shifting of the fast system energy-levels 
by a combination of Berry phase and stochastic effects; and (2) the semiclassical limit of 
the slow system effective dynamics obeys a Langevin equation in which geometrical effects 
produce Lorentz-like and electric-like forces. 

In the semiclassical and low- noise limit, we calculate approximately the average level- 
shift and broadening produced by this geometric mechanism. Monte Carlo evaluation of 
eqns. (|50|) and (|5T|) would be very interesting as this would free the analysis from the low- 



noise limit. 

Formally, the semiclassical limit was taken by letting Tl — > 0. In fact, h is finite, and 

an experimental realization of this limit must be approached differently. Formally, the 

semiclassical limit is controlled by the real exponential factor appearing in eqn. (p8|), and 

is approached when it becomes sharply peaked. This occurs when h — > 0, though more 

generally when r/R kTT ^> h 2 . Thus strong damping is one way to produce semiclassical 

2 

behavior. The low- noise limit required kT <C r/R /T. Thus, sufficiently low temperature 
will produce low-noise behavior. Our approximate results for the average level-shift and 
broadening are thus expected to apply in the limit of strong damping and sufficiently low 
temperature. 
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In an interesting Comment, Simon and Kumar |2T| propose a physical setting in which 
Berry's phase should produce level-broadening. Their underlying idea is similar to the one 
we propose: a range of level-shifts occur producing a broadening of the original energy level. 
They do not provide a formal development of their proposal, however. Gamliel and Reed 
|22|| consider the original Berry phase scenario of a pseudo-spin (fast system) interacting 
with an adiabatically evolving pseudo- magnetic field (slow system). However, they allow 
for the presence of a stochastic process whose sole effect is to relax the pseudo-spin to an 
equilibrium state (which might evolve with time). The stochastic process is assumed to 
produce no Berry phase effects in the pseudo-spin dynamics. The deterministic motion of 
the pseudo-magnetic field produces a unique level-shift in each fast system energy 

level E. The central question for these authors is whether the shift 5i?[R] is observable in 
the presence of the conventional level-broadening produced by the relaxation process. In 
our scenario, the slow system motion is the random process, and the level-broadening arises 
through the Berry phase induced level-shifts. 
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APPENDIX: 

In this Appendix we consider a free classical Brownian particle moving in one-spatial 
dimension. We will establish a condition that insures that the Brownian motion is adiabatic 
relative to a quantum system that interacts with the Brownian particle. 

The Brownian motion is described by Langevin's equation, 

mR + r)R = F(t), (Al) 

where m is the particle mass; rj is the friction coefficient; and F(t) is a stochastic force with 
correlation function, 

( F(t)F(f) ) = 2rjkT5(t - t') . (A2) 

The stochastic force F(t) fluctuates rapidly. The duration r c of a force fluctuation is known 
as the correlation time. It is determined by the time scale of the microscopic processes which 
produce the stochastic force. The particle velocity R varies on a much slower time scale due 
to the large inertia of the Brownian particle. 

Since the Brownian particle is free, E = mv 2 /2. The time-averaged energy is 

(E) = -(v 2 ) = -KM, 

where K v (t) = ( v(0)v(t) ) is the velocity correlation function. Its Fourier transform is the 
spectral density J v (lu) of v(t) f23fl , 
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K v (t) = / diO COSLUt J v {lo) . 
JO 



Thus, 



m r°° 

(E) = - dujJ v (u) . (A3) 
2 jo 

Defining the energy spectral density Pe{^) as the mean noise-energy in the frequency range 
(a;, oj + du), we see from eqn. ( |A3| ) that, 

in 

p E {u) = -J v {u) . (A4) 

It is well-known [^3| that, J v {uS) = limr-+oo 2|v(tf)| 2 /T. Since v(u) = —iojR(uj), eqn. (|A4|) 
becomes, 

Pe\w) = —^-J r (oj) , (A5) 

where J R (u) = \im T ^ 00 2\R(uj)\ 2 /T. 

It is preferable to express Pe{^->) in terms of Jp(u). To do this we use the Langevin 
equation to relate Jr{uj) to Jp(ou). Fourier transforming eqn. (|A1|) gives, 

-F{u) 



R(lu) 



where u r = r]/m. Thus, 



and, 



muj (uj + iuj r 



m 2 uj 2 (oj 2 + uj 2 ) 



2m (oj 2 + u 2 ) 



In writing the force correlation function in eqn. ( |A2|) as being proportional to 5(t — t'), we 
assumed r c was effectively zero, and Jp(u) = 2r]kT for all oj. In fact, r c is not zero, though 
macroscopically small. Thus Jf(oj) = 2rjkT only up to a cutoff frequency A ~ l/r c so that 

= /fl 6(A - W ) , (A6) 
m (ur + uj;) 



where Q(x) vanishes if x < and is 1 otherwise. From eqns. ( |A3[ ) — ( |A5| ) , and ( |A6[ ) we see 
that noise fluctuations in the random process R(t) only exist for frequencies oj < l/r c . If the 
Brownian particle is coupled to a quantum system whose energy-level spacing is AE, the 
Brownian motion will be unable to produce transitions in the quantum system if fi/r c <C AE. 
If this condition is satisfied, the quantum system will see the Brownian motion as adiabatic. 
This is the desired adiabaticity condition. 
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